Measurement of Direct Photons in Au+Au Collisions at ^/s^ = 200 GeV 
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Wo report the measurement of direct photons at midrapidity in Au+Au collisions at ys„jy — 200 
GcV. The direct photon signal was extracted for the transverse momentum range of 4 GeV/c < pt < 
22 GeV/c, using a statistical method to subtract decay photons from the inclusive photon sample. 
The direct-photon nuclear-modification factor Raa was calculated as a function of pr for different 
Au+Au collision centralities using the measured p+p direct-photon spectrum and compared to 
theoretical predictions. Raa was found to be consistent with unity for all centralities over the 
entire measured pr range. Theoretical models that account for modifications of initial-direct-photon 
production due to modified-parton-distribution functions (PDFs) in Au and the different isospin 
composition of the nuclei predict a modest change of Raa from unity. They are consistent with 
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the data. Models with compensating effects of the quark-gluon plasma on high-energy photons, 
such as suppression of jet-fragmentation photons and induced-photon bremsstrahlung from partons 
traversing the medium, are also consistent with this measurement. 

PACS numbers: 25.75.Dw 



Direct photons are a powerful probe to study ultra- 
relativistic heavy-ion collisions where a hot and dense 
quark-gluon plasma (QGP) is formed. Direct photons 
are defined as all photons that arise from processes during 
the collision, rather than from decays of hadrons. The 
biggest challenge in the measurement of direct photons is 
to distinguish them from the large background of decay 
photons. 

Direct photons with intermediate and high transverse 
momentum (px > 4 GeV/c) are produced predomi- 
nantly from initial-hard-scattering processes of the collid- 
ing quarks or gluons, such as q+g q+j or q+q .9+7- 
In addition, they can be produced as bremsstrahlung 
emitted by a scattered parton, from the fragmentation 
of quarks and gluons, or from the interaction of a scat- 
tered parton with the medium created in heavy-ion col- 
lisions [ll-ll]- Additional photons may be emitted at low 
transverse momentum as thermal radiation from the par- 
tonic and hadronic phases. 

The production of direct photons, if compared to the 
scaled rates, is also affected by possible modifications 
of the initial state of the colliding nuclei, like shadowing 
and ant i- shadowing, and by the different isospin compo- 
sition of Au nuclei in contrast to protons, as explained 
in Sec. 2.2 of 0. In addition, the different quark charge 
squared content of p and n influences the yields from ini- 
tial hard scattering in heavy ions, as explained in detail 
in Sec. 3.3 of Q. 

The direct photons should not be affected by the 
medium as they traverse it, since they are both electri- 
cally and color neutral, but the presence of the medium 
can affect the total direct photon yield. For instance, 
parton energy loss 0-11] can reduce the fraction of frag- 
mentation photons at a given pT, while the scattering 
of a hard parton on a thermal one can produce a high 
Pt photon with approximately the same momentum as 
the original parton (jet-photon conversion) [if. As a re- 
sult, theoretical models predict that the yield of direct 
photons in Au-|-Au collisions will be somewhat modified 
compared to the scaled yield from p+p collisions at the 
same energy [l-IHj. 

Previous PHENIX measurements of direct photon 
spectra in Au-|-Au at = 200 GeV, from the 2002 

RHIC run, showed no significant deviation above pt > 
6GeV/c from the scaled invariant yield of NLO pQCD 
predictions for p+p collisions [l^l- On the other hand, 
PHENIX measurements of virtual direct photons at low 
transverse momentum {px < 4GeV/c) found a large ex- 
ponentially distributed excess of direct photons, com- 
pared to the scaled p+p QCD prediction, which was 



attributed to thermal photon radiation from the QGP 
formed in central Au-|-Au collisions 11]. Although the 
contribution of thermal photons at large px quickly di- 
minishes due to the exponential falloff, the direct photon 
yield at high px may be modified in nuclear collisions 
due to the various effects mentioned above and deserves 
careful study. 

We report on the measurement of direct photons in 
Au+Au colhsions at y/s^ = 200 GeV at RHIC from 



data taken by the PHENIX experiment [IJ] in 2004. The 
analysis used 1.03 x 10^ minimum bias events, which is 
more than a tenfold increase compared to the previous 
measurement [lo| . The centrality (impact parameter) of 
the Au-|-Au collision was determined from the correla- 
tion between the number of charged particles detected in 
the Beam-Beam Counters (BBC), in the pseudorapidity 
range 3.0 < \ri\ < 3.9, and the energy measured in the 
zero-degree calorimeters (ZDC). The average number of 
binary nucleon-nucleon collisions {{Ncoii}) was estimated 
for each centrality bin with a Glauber Model Monte Carlo 
that simulated the BBC and ZDC responses 13]. 

Photon candidates were reconstructed in the electro- 
magnetic calorimeter (EMCal) located in the central 



arms of PHENIX [iJ. The EMCal covers \ri\ < 0.35. 



It comprises six sectors of lead-scintillator calorimeter 
(PbSc) and two sectors of lead-glass Cerenkov calorime- 
ter (PbGl). Located at a radial distance of about 5 m, 
the two subsystems cover a total of tt in azimuth. The 
segmentation is x Ary ~ 0.011 x 0.011 for PbSc and 
- 0.008 X 0.008 for PbGl. 

At high transverse momenta, the minimum opening 
angle between the two photons of a 7r° decay decreases, 
and the distance between the two clusters on the EMCal 
surface becomes comparable to the tower segmentation, 
with the result that the showers begin to merge. This 
starts to occur at ^lOGeV/c (16GeV/c), and affects 
50% of all tt" decays at px -16GeV/c (24GeV/c) for 
the PbSc (PbGl) detector. Furthermore, when the decay 
photons partially overlap, but are not yet indistinguish- 
ably merged, the energy may be imperfectly shared by 
the clustering algorithm: one reconstructed cluster has 
more, the other has less energy than the original photons. 
This may change the apparent photon yield, particularly 
at high Px- The effect is significant for the PbSc, but 
negligible for the PbGl in the measured px range. 

Direct photons were measured on a statistical basis, as 
in earlier PHENIX direct photon measurements [13, [3] ■ 
Photon-like clusters were identified by applying Parti- 
cle Identification (PID) cuts based on the parameterized 
shower profile for a photon. The analyses were performed 



4 



independently with the PbSc and the PbGl calorimeters, 
and the fully corrected results were combined. Since the 
methods were different for the two detectors, the system- 
atic uncertainties are uncorrelated. 

In the analysis of the PbGl data, the w 10 — 15 % con- 
tamination of the photon candidate spectra with charged 
particles was subtracted by associating photon candi- 
dates with charged hits in the pad chamber (PCS) sit- 
uated directly in front of the calorimeter. Neutrons 
and anti-neutrons (5% contribution to the cluster energy 
spectrum at 4 GeV, vanishing at 7 GeV) were subtracted 
based on a full GEANT [l^ simulation of the detector 
response. 

The spectra were also corrected for the loss of photons 
due to conversions into e~^e~ pairs in the material in front 
of the PbGl. The resulting spectra were corrected for the 
detector acceptance and reconstruction efficiency. The 
acceptance is influenced by the detector geometry and 
the exclusion of detector areas from the analysis. It was 
calculated with a Monte Carlo simulation. The efficiency 
correction takes into account the energy resolution of the 
calorimeter, the applied PID cuts, and occupancy effects 
in the high- multiplicity environment. The reconstruc- 
tion efficiency was determined by embedding simulated 
photons into real events and analyzing these embedded 
photons with the same analysis cuts. 

Merged clusters at high pT were removed by the PID 
cuts. Photons from hadron decays, mainly from 7r° and rj, 
measured by PHENIX [l3, [3| , were simulated using the 
decay kinematics and detector geometry and, following 
the method used in earlier analyses [l^l , used to calculate 

data j data 

the ratio = '''"■jfer"; o = ^"i^l^f'"" . This ratio 

^decayl^MC ^ decay 

was used to extract the direct photon invariant yield via 

^direct — (1 "j^^^inclusive ■ 

In the PbSc analysis, photon candidates (clusters pass- 
ing PID cuts) were corrected for the fraction of electrons, 
charged hadrons, and neutrons passing those PID cuts; 
this fraction was derived from full GEANT detector simu- 
lations using particle spectra measured by PHENIX. The 
result was the raw inclusive photon distribution. Note 
that at higher pr the calorimeter response to true single 
photons and correlated decay photons is different, there- 
fore, the raw inclusive spectrum cannot be corrected in 
one simple step. Instead, first the expected raw distri- 
bution of background photons from hadron decays (pre- 
dominantly tt" and 77), containing all detector effects, de- 
noted '^'Jdicay.raw^ Calculated in a full GEANT simula- 
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tion which used the measured tt*^ and 77 spectra 
After subtracting the raw decay photons from the raw in- 
clusive photons, the acceptance and efficiency correction 
for the remaining direct (single) photons was obtained 
using simulated single photons embedded in real events. 
These corrections were then applied to the raw direct 
single-photon distribution to get the final direct-single- 
photon distribution J^f^^^f, the experimental result. The 



TABLE I: Statistical and systematic uncertainties of the di- 
rect photon yield, in %, for Au-f An minimum bias events, as 
estimated for the measurement with the PbGl (PbSc). The 
values for three transverse momenta are given. All systematic 
uncertainties are correlated in pr- 



Error type / pr 


4.75 GeV/c 


9.25 GeV/c 


15GeV/c 


Background corrections 


9.1(5.2) 


5.7(2.5) 


5.1(2.2) 


Yield corrections 


11.9(10.5) 


8.3(9.4) 


7.9(11.2) 


Energy scale 


7.9(6.8) 


6.8(7.0) 


6.8(7.0) 


Decay 7 simulation 


12.5(7.2) 


5.2(4.3) 


3.8(3.7) 


Total Systematic 


21.0(13.9) 


13.2(12.7) 


12.3(13.9) 


Total Statistical 


0.9(0.4) 


4.1 (2.6) 


8.8(8.2) 


Combined Systematic 


11.6 


9.1 


9.3 


Combined Statistical 


0.4 


2.1 


5.9 



result from the 27^^^ in the PbSc analysis cannot 
be simply acceptance corrected to get the true decay-7 
background, therefore, the ratio for the PbSc was cal- 
culated using the decay photon background Monte Carlo 

^data +^MC 

calculation from the PbGl analysis as Rj = "^"""mc ■ 

For both the PbGl and PbSc measurements, there are 
four distinct sources of systematic uncertainties (TableHl, 
all of which are pr-correlated. Uncertainties from back- 
ground corrections come from the subtraction of hadron 
and electron contamination and corrections for photon 
conversions. The corrections to the raw yields by the sim- 
ulations are another source of uncertainty. The energy 
scale of the calorimeters is only known with a 1.2% preci- 
sion and thus leads to uncertainties in the direct photon 
measurement. The decay photon calculation adds fur- 
ther systematic uncertainties due to the extraction of the 
.0 a the parameterization of the input hadron spectra 
and ratios such as t^/tt". 

The fully corrected results obtained with the two anal- 
yses agree within their respective uncertainties and were 
combined. The spectra and nuclear modification fac- 
tors in this publication represent a weighted average of 
the two independent measurements. Since the system- 
atic uncertainties are taken to be uncorrelated between 
the two analyses, the weight w was determined from 
their total uncertainty aTotai, which is the quadratic 
sum of all statistical and systematic uncertainties, by 

The ratio R-y is shown in Fig. [1] for minimum bias 
Au-|-Au collisions and the two extreme centrality bins. 
An excess above unity indicates the presence of direct 
photons. Such an excess is clearly visible for all central- 
ity selections. The ratio R^ increases with centrality due 
to the suppression of the 7r° 0, |1] and the associated 
decay photons. 

The combined direct photon spectra in Au-f-Au colli- 
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FIG. 1: Ratio R-y for difFerent centrality selections, for the 
PbGl and the PbSc analysis. The error bars indicate point- 
to-point uncertainties, the boxes around the points indicate 
Pt correlated uncertainties. 
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sions are shown in the top panel of Fig.[2]for ten centrality 
selections. The shape of the spectra are seen to be similar 
for all centralities. The bottom panel shows a compari- 
son of the PbGl and PbSc spectra to the combined result 
for the — 5% most central collisions. A good agreement 
between the two measurements is observed. 

Fig. [5] also includes the p+p spectrum at the same en- 
ergy, measured by PHENIX [20|. The p+p spectrum 
is compared to a power law fit {A/px)" with power 
n = 7.08 ± 0.09(stat) ± O.l(syst) obtained by fitting the 
region pr > S GeV/c (l^l- The fit is extrapolated to 
lower Pt- A power law fit to the minimum bias (most 
central) Au-|-Au spectrum yields a power of n = 6.85 ± 
0.07(stat)±0.02(syst) {n = 7.18±0.14(stat)±0.06(syst)) 
consistent with the power of the p+p fit. The agreement 
indicates no apparent shape modification of the spectra 
compared to p+p collisions. 

For hard processes, the yield in A-l-A collisions for a 
particular impact parameter selection is expected to be 
equal to the cross section in p+p collisions, scaled by the 
average nuclear thickness function (Taa) = (-^coii) /fpp'^' 
for the associated centrality selection. Here, (Ncow) is the 
number of binary nucleon-nucleon collisions, calculated 
with the Glauber Model Monte Carlo for the selected 
centrality, and (J^^^ is the total inelastic p+p cross sec- 
tion of 42 mb. In Fig. [51 the power law fit to the p+p 
direct photon spectrum has been scaled by the nuclear 
thickness function for each of the ten centrality selections, 
and overlaid on the measured result for that centrality. 
The comparison indicates that the magnitude, as well as 
the shape of the direct photon spectra, are in agreement 
with expectations from p+p collisions for all centralities. 

Nuclear effects are quantified by the nuclear modifica- 
tion factor, i?AA- For a given centrality selection, Raa 
is given by the ratio of the measured invariant yields in 



FIG. 2: (a) Direct photon spectra for all centrality selections 
in Au-|-Au, and for p+p measured in [S^]. The error bars 
indicate point-to-point uncertainties, the boxes around the 
points indicate pr correlated uncertainties. The lines depict 
a Taa scaled fit for pr > 8GeV/c to the p+p cross section, 
they are dashed for the range where the fit is extrapolated to 
lower Pt- (b) Comparison of the PbGl and PbSc results with 
the combined result for the — 5% most central events. The 
error bars show the total uncertainties. 



Au-|-Au collisions, divided by the production cross sec- 
tion for the same particle in p+p collisions, scaled with 
the average nuclear thickness function for that centrality: 



Raa{pt) 



{l/Nfl)<fNAAldpTdy 
(Taa) X (Papp/dprdy 



(1) 



where (Pupp/ dprdy is the measured p+p cross section for 
direct photons (20j . 

The direct photon nuclear modification factor is shown 
in Fig.[3]for three different centrality selections. The i?AA 
results are calculated using the measured direct photon 
results from p+p collisions for the first time. The -Raa 
values are consistent with unity, within errors, for all 
centrality selections over the entire px range. 

In Fig. m the measured nuclear modification factor for 
central Au-|-Au collisions is compared to theoretical cal- 
culations that predict modifications of the direct pho- 
ton yield due to initial state (IS) and final state (FS) 
effects 043- IS effects include the isospin effect due to 
the difFerent photon cross sections in p+p, n+n, and p+n 
collisions ("Isospin effect" in Fig. 0]), and modifications 
of nuclear structure functions due to shadowing and anti- 
shadowing ("EPS09 PDF") Q. The EPS09 calculation 
also includes the isospin effect. 

FS modifications due to QGP lead, on one hand, 
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to a lower photon yield, since energy loss of a parton 
also means suppression of the corresponding fragmen- 
tation photon yield. On the other hand, QGP effects 
can increase the photon yield due to radiation resulting 
from jet- medium interactions ( "prompt -l-QGP" ) 0, 
This FS calculation also takes into account the afore- 
mentioned IS effects. Yet another calculation Q in- 
cludes IS effects, as well as FS energy loss and medium- 
induced photon bremsstrahlung and the LPM effect 
( "coherent -|-conversion-|- AE" ) . The data are consistent 
with a scenario where the hard scattered photons are 
produced taking account of the isospin effect and mod- 
ifications of the nuclear PDFs and then simply traverse 
the matter unaffected. Balancing effects from the QGP 
such as fragmentation photon suppression and enhance- 
ment due to jet-medium interactions are not excluded by 



the data. The approach in [sj is in disagreement with 
the data. This confirms that the majority (if not all) 
direct photons at high px come directly from hard scat- 
tering processes and suggests that possible effects from 
the QGP all but cancel. 

In summary, PHENIX has measured direct photon 
spectra in Au+Au coUisions at y/s^ = 200 GeV at 
midrapidity in the transverse momentum range of 4 < 
Pt < 20GeV/c. The direct photon nuclear modification 
factor Raa has been calculated as a function of pT us- 
ing a measured p+p reference for the first time. It is 
consistent with unity for all centrality selections over the 
entire measured px range. Theoretical models for direct 
photon production in Au-I-Au collisions are compared to 
the data. Some of these models are found to be in quan- 
titative agreement with the measurement while others 
appear to be disfavored by the data. Collectively, the ef- 
fects of the QGP on the high pT direct photon yield are 
apparently small. 

We thank the staff of the Collider- Accelerator and 
Physics Departments at Brookhaven National Labora- 
tory and the staff of the other PHENIX participating in- 
stitutions for their vital contributions. We acknowledge 
support from the Office of Nuclear Physics in the Office 
of Science of the Department of Energy, the National Sci- 
ence Foundation, Abilene Christian University Research 
Council, Research Foundation of SUNY, and Dean of 
the College of Arts and Sciences, Vanderbilt University 
(U.S. A), Ministry of Education, Culture, Sports, Science, 
and Technology and the Japan Society for the Promotion 
of Science (Japan), Conselho Nacional de Desenvolvi- 
mcnto Cientifico e Tecnologico and Fundagao de Amparo 
a Pesquisa do Estado de Sao Paulo (Brazil), Natural Sci- 
ence Foundation of China (P. R. China) , Centre National 
de la Recherche Scientifique, Commissariat a I'lilnergie 
Atomique, and Institut National de Physique Nucleaire 
et de Physique des Particules (France), Ministry of In- 
dustry, Science and Tekhnologies, Bundesministerium fiir 
Bildung und Forschung, Deutscher Akademischer Aus- 
tausch Dienst, and Alexander von Humboldt Stiftung 
(Germany), Hungarian National Science Fund, OTKA 
(Hungary), Department of Atomic Energy (India), Israel 
Science Foundation (Israel), National Research Founda- 
tion and WCU program of the Ministry Education Sci- 
ence and Technology (Korea) , Ministry of Education and 
Science, Russian Academy of Sciences, Federal Agency of 
Atomic Energy (Russia), VR and the Wallenberg Foun- 
dation (Sweden), the U.S. Civilian Research and Devel- 
opment Foundation for the Independent States of the 
Former Soviet Union, the US-Hungarian NSF-OTKA- 
MTA, and the US-Israel Binational Science Foundation. 



Deceased 



7 



' PHENIX Spokesperson: jacak@skipper.physics.sunysb.edu 
[1] R. J. Fries, B. Muller, and D. K. Srivastava, Phys. Rev. 

Lett. 90, 132301 (2003). 
[2] S. Turbide, C. Gale, E. Frodermann, and U. Heinz, Phys. 

Rev. C 77, 024909 (2008). 
[3] I. Vitev and B.-W. Zhang, Phys. Lett. B 669, 337 (2008). 
[4] C. Gale, arXiv:0904.2184 (2009). 

[5] F. Arleo, K. J. Eskola, H. Paukkunen, and C. A. Salgado, 

JHEP 1104, 055 (2011). 
[6] F. Arleo, JHEP 09, 015 (2006). 

[7] K. Adcox et al. (PHENIX Collaboration), Phys. Rev. 

Lett. 88, 022301 (2002). 
[8] S. Adler et al. (PHENIX Collaboration), Phys. Rev. Lett. 

91, 072301 (2003). 
[9] X.-N. Wang and M. Gyulassy, Phys. Rev. Lett. 68, 1480 

(1992). 

[10] S. S. Adler et al. (PHENIX Collaboration), Phys. Rev. 

Lett. 94, 232301 (2005). 
[11] A. Adare et al. (PHENIX Collaboration), Phys. Rev. 

Lett. 104, 132301 (2010). 



[12] K. Adcox et al. (PHENIX Collaboration), Nucl. Instrum. 

Meth. A 499, 469 (2003). 
[13] M. L. Miller, K. Reygers, S. J. Sanders, and P. Steinberg, 

Annu. Rev. Nucl. Part. Sci. 57, 205 (2007). 
[14] L. Aphecetche et al. (PHENIX Collaboration), Nucl. In- 
strum. Meth. A 499, 521 (2003). 
[15] S. S. Adler et al. (PHENIX Collaboration), Phys. Rev. 

Lett. 98, 012002 (2007). 
[16] GEANT 3.2.1, CERN Program Library (1993), 

|http: //wwwasdoc .web. cern. ch/wwwasdoc/pdf dir/geant -pdf] 
[17] A. Adare et al. (PHENIX Collaboration), Phys. Rev. 

Lett. 101, 232301 (2008). 
[18] A. Adare et al. (PHENIX Collaboration), Phys.Rev. 

C82, 011902 (2010). 
[19] K. Nakamura et al. (Particle Data Group), J. Phys. G 

37, 075021 (2010). 
[20] A. Adare et al. (PHENIX Collaboration), 

arXiv: 1205.5533 (2012), submitted to Phys. Rev. 

D. 



